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Abstract 



We present the complete NLO QCD analysis of the photon structure 
functions F^(x, Q 2 ) and F2(x, Q 2 ) for a real photon target. In particular we 
study the heavy flavor content of the structure functions which is due to two 
different production mechanisms, namely collisions of a virtual photon with 
a real photon, and with a parton. We observe that the charm contributions 
are noticeable for F 2 7 (x, Q 2 ) as well as F2(x, Q 2 ) in the x-region studied. 
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1 Introduction 



In the past two decades there has been considerable interest in the study of 
photon-photon interactions in electron-positron colliders. When one photon 
is virtual and the other one is almost real the analogy with deep-inelastic 
electron-nucleon scattering motivates the introduction of the corresponding 
structure functions F^(x,Q 2 ) (k = 2,L) for the photon. The deep-inelastic 
structure function F2(x, Q 2 ) was originally measured by the PLUTO col- 
laboration H at PETRA using single-tag events in the reaction e~ + e + — > 
e~ + e + + hadrons. In the past several years there has been a series of new 
measurements at PETRA, PEP and TRISTAN by several groups, includ- 
ing CELLO 0, TPC2 7 0, TASSO |J, JADE §, AMY §, VENUS 
and TOPAZ ||. All these groups concentrated on the measurement of the 
light-quark contribution to F^x^Q 2 ). The heavy-quark component (mainly 
charm) has been hard to extract due to problems identifying charmed parti- 
cle decays so its contribution to the data was sometimes removed according 
to a Monte Carlo estimate. In the near future higher-luminosity runs at 
TRISTAN should yield some information on heavy-quark (mainly charm) 
production and this is one reason that we study it here. At this moment 
the available data for F^(x, Q 2 ) are in the region 0.03 < x < 0.8 and 
1.31 (GeV/c) 2 < Q 2 < 390 (GeV/c) 2 . Due to the experimental limitation 
that xy 2 « 1 (for a definition of x and y see (2.5)), there are no data avail- 
able for the longitudinal structure function Fl(x,Q 2 ). However there exists 
some hope that F2(x, Q 2 ) can be measured |§ at LEP. Finally two-photon 
reactions are important to understand as background processes to the normal 
s-channel reactions at present and future e + e~ colliders. These machines will 



have a large amount of beamstrahlung |TIJ , [ TT| . Therefore a basic input is 
the parton density in a photon which will be modified if higher order pQCD 
corrections are included. 

As far as theory is concerned the first attempt to give a theoretical de- 
scription of the photon structure function in the context of perturbative 



QCD was given by E. Witten in |fI2| . He suggested that both the x and 
the Q 2 dependence of these structure functions were calculable in pQCD at 
asymptotically large Q 2 . Thus from a theoretical point of view this process 
should provide a much more thorough test of pQCD than the corresponding 
deep- inelastic scattering off a nucleon target, where only the Q 2 evolution of 
the structure functions is calculable. The original optimism subsided once it 
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was realized that there were complications with experimental confirmation of 
this prediction at experimental (non- asymptotic) values of Q 2 [13] , ||14j| . For 



recent reviews see |T5[] . In particular there is a contamination of the purely 
pointlike pQCD contribution by the hadronic component of the photon. This 
latter piece, which is most important at small virtualities, is not calculable 
in pQCD and must be extracted from experimental data. One of the ap- 
proaches used is to describe this hadronic piece by parton densities in the 
photon, analogous to the parton densities in a hadronic target. For parame- 
terizations see il6f , fl7| , [jiS|] , |19| , p0| and plfl . For a different approach see 

In |19| a next to leading order (NLO) analysis was carried out for the 
photon structure function F^ix, Q 2 ). This analysis also includes the lowest 
order contribution coming from heavy flavor production, which is described 
by the Bethe-Heitler cross section corresponding to the process 7* + 7 — > 
Q + Q. In this case the mass m of the heavy flavor is not neglected with 
respect to Q 2 especially in the threshold region. If Q 2 >> m 2 one encounters 
large logarithmic terms containing ln(Q 2 /m 2 ), which have to be summed 
using the Altarelli-Parisi (AP) equations. This procedure provides us with 
the heavy flavor densities in the photon which are akin to the parton densities 
originating from the light quarks in the photon. The same procedure has been 
applied for the longitudinal structure functions F2(x, Q 2 ) in [Q but only in 
leading order. 

In this paper we want to extend the above analysis by including higher 
order pQCD corrections which were not considered in the literature so far. 
Since the NLO QCD corrections to the longitudinal coefficient functions due 
to massless partons |24| and heavy flavors have been recently calculated 
we are now also able to present a NLO analysis for F2(x,Q 2 ). In addition 
we can also improve our knowledge of the heavy flavor content of F^(x, Q 2 ) 
by including the order a s corrections to the Bethe-Heitler process 7* + 7 — ► 
Q+Q. We also include corrections to F^(x, Q 2 ) (k = 2, L) due to heavy flavor 
production mechanisms given by the processes 7* + g — > Q + Q (corrected 
up to order a 2 ) and 7* + q(q) — > Q + Q + q(q), where the incoming gluon 
and (anti) quark originate from the on- mass-shell photon. Furthermore we 
use the most recent gluon and (anti) quark densities in our analysis. 

Finally we should mention that there was a previous investigation of 
pQCD corrections to heavy quark production in [26], where it was assumed 



that both photons were off-mass-shell and a small value for the photon vir- 
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tuality was chosen for generating numerical results. Since these authors did 
not therefore encounter mass singularities they had no need to perform any 
mass factorization. Hence their method was different from the one we adopt. 

The paper is organized as follows. In section 2 we present the photonic 
and hadronic coefficient functions corrected up to next to leading order in 
a s , which are needed for the photon structure functions F^(x, Q 2 ) (k = 2,L). 
In section 3 we show the differences between the leading order (LO) and the 
next to leading order (NLO) photon structure functions. In particular we 
discuss the effect of the heavy flavor component (mainly charm) originating 
from the hadronic as well as the pointlike photon interactions. 
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2 Higher-Order Corrections to the Photon 
Structure Functions 



The deep-inelastic photon structure functions denoted by F^(x,Q 2 ) {k = 
2,L) are measured in e~e + collisions via the process (see fig.l) 

e~(p e ) + e + ^e-(p' e )+e + + X , (2.1) 

where X denotes any hadronic state which is allowed by quantum-number 
conservation laws. When the outgoing electron is tagged then the above 
reaction is dominated by the photon-photon collision reaction (see fig.l) 

7*{q)+7(k)^X, (2.2) 

where one of the photons is highly virtual and the other one is almost on- 
mass-shell. The process (2.1) is described by the cross section 

d 2 a f , S x 27ca 2 S 



J dzzf*(z 



dxdy J ' r 'm^ Q 4 

{1 + (1 - y) 2 }F?(x, Q 2 ) - y 2 F2(x, Q 2 )} , (2.3) 



x 



where (x, Q 2 ) (k = 2, L) denote the deep-inelastic photon structure func- 
tions and a = e 2 /4ir is the fine structure constant. Furthermore the off-mass- 
shell photon and the on-mass-shell photon are indicated by the four-momenta 
q and k respectively with q 2 = —Q 2 < and k 2 ~ 0. Because the photon 
with momentum k is almost on-mass-shell, expression (2.3) is written in 
the Weizsacker- Williams approximation. In this approximation the function 
f^(z, S/m 2 ) is the probability of finding a photon ■y(k) in the positron, (see 
fig.l). The fraction of the energy of the positron carried off by the photon is 
denoted by z while v^S* is the cm. energy of the electron-positron system. 
The function f^(z, S/m 2 ) is given by (see f27j) 

/ e fe A ) ^i±(^£)! ln a-)(f 2 -^) , (2 . 4) 

provided a heavy quark with mass m is produced. The scaling variables x 
and y are defined by 

Q 2 k-q , 

x = — 1 — > y = -, — > q = Pe-p e , (2.5) 

2k ■ q k ■ p e 
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where p e , p' e are the momenta of the incoming and outgoing electron respec- 
tively. Following the procedure in [28 the photon structure functions in the 
QCD-improved parton model have the following form 




where the meaning of the symbols is explained below. 

The quantities £ 7 and A 7 represent the singlet and non-singlet combina- 
tions of the quark densities in the photon respectively while the gluon density 
is represented by g 1 '. The same flavor decomposition is also applied to the 
hadronic (Wilson) coefficient functions Ck,% (i = q,g) so that C^ q (z,Q 2 /M 2 ) 
and C^ q s (z,Q 2 /M 2 ) stand for the singlet and non-singlet coefficient func- 
tions respectively, and Ck, g (z,Q 2 /M 2 ) denotes the gluonic coefficient func- 
tion, where M 2 is the mass factorization scale. The hadronic coefficient 
functions can be attributed to hard processes with a light quark or gluon 
in the initial state, such as 7* + q — > q + g or 7* + g — > q + q, where the 
initial parton emerges from the real (on-mass-shell) photon. Hence they are 
multiplied by the corresponding parton densities in the photon. 

We also make a distinction between light and heavy flavor contributions 
to the coefficient functions. The latter are indicated by their explicit de- 
pendence on the heavy flavor mass m. For example in the contribution to 
Ck,i(z, Q 2 /M 2 , m 2 ) (second part of (2.6)) the virtual photon is attached either 
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to the incoming light quark as is the case in the reaction , y* + q^q + Q + Q 
or indirectly to the incoming gluon. Actually the Ck t i(z,Q 2 /M 2 ,m 2 ) be- 
long to the same class as the hadronic light parton coefficient functions 
presented in the first part of expression (2.6). The only difference is that 
Ck,i( z ,Q 2 /M 2 ,m 2 ) receives contributions from a heavy flavor pair produced 
in the final state. 

In the third set of terms in (2.6) the heavy flavor coefficient functions 
originate from subprocesses where the virtual photon is attached to one of the 
outgoing heavy flavors, as for example in 7*+g — > Q + Q, so they are given an 
additional superscript H. Finally the fourth set of terms in (2.6) contain the 
photonic coefficient functions indicated by Cfc i7 coming from reactions such 
as7* + 7^g + gor7*+7-H><5 + <5- These originate from hard processes 
where the (on-shell) real photon is directly attached to the light or heavy 
quarks produced in the final state so there is no need for any convolution 
integral. 

The index % in (2.6) runs over all light active flavors whose number is given 
by 71/ and e,, e# stand for the charges of the light and heavy quarks respec- 
tively in units of e. The upper boundary of the integrals in (2.6) containing 
the convolution of the parton densities with the heavy flavor coefficient func- 
tions is given by 

Zmax = 4m 2 + Q 2 ■ (2,7) 

The parton densities as well as the coefficient functions depend on the mass 
factorization scale M except for the Cj?^ which can be calculated in pQCD 
without performing mass factorization. Notice that in addition to the mass 
factorization scale M the quantities in (2.6) also depend on the renormaliza- 
tion scale R which appears in the pQCD corrections via a s (R 2 ). However in 
this paper we will put R = M. 

According to the origin of the photonic parton densities and the two differ- 
ent types of coefficient functions i.e., Ck, q , Ck, g (hadronic) and Ck a (photonic) 
we will call the first three terms in (2.6) (represented by the integrals), the 
hadronic photon parts, and the last term the pointlike photon part. Notice 
that both these terms are separately factorization scheme dependent as indi- 
cated by the presence of the scale M. In particular the scheme dependence 
of the pointlike photon part in (2.6) is due to the light quark contribution 
Ck, 7 (x, Q 2 /M 2 ). The scheme dependence is cancelled by the hadronic photon 



8 



part due to the light quark contribution provided that the quark densities 
and the hadronic coefficient functions are computed in the same scheme as 
Cfc i7 (x, Q 2 /M 2 ). The hadronic heavy flavor part is scheme dependent in it- 
self. The photonic heavy flavor piece is obtained without having to perform 
mass factorizaton and needs no parton distribution functions, and is thus not 
dependent on the factorization scheme. 

In the subsequent part of this section we will discuss the contributions 
to the coefficient functions in (2.6) which are needed for a next to leading 
order (NLO) description of the photon structure functions Fl (x, Q 2 ) and 
F2(x,Q 2 ). The results of our calculations will be presented in the plots of 
section 3. For these NLO calculations we also have to use the next to leading 
logarithmic (NLL) approximation to the parton densities, which are given 
for example in p|,[|lj,|21|. 



Starting with the NLL parton densities the singlet and nonsinglet combi- 
nations are written in the following way. Below the charm-quark threshold 
we have 

»/ = 3 , XX 2 = | , (2-8) 
•=i 6 i=i 9 

= v? + u 1 + <F + dT + s 7 + s 7 , (2.9) 
A 7 = i(2w 7 + 2u< - d? - J 7 - s 7 - s 1 ) . (2.10) 

Above the charm-quark threshold and below the bottom-quark threshold the 
above quantities are changed into 

4 in 4 34 

i=i 9 i=i 81 
S 7 = u< + u 1 + (T + J 7 + s 7 + s 7 + c 7 + c 7 , (2.12) 

A 7 = i(w 7 + m 7 + c 7 + c 7 - d? - J 7 - s 7 - s 7 ) . (2.13) 
6 

Finally above the bottom-quark threshold they become 

n/ = 5 , Ee, 2 = ^ , E«? = i' ( 2 - 14 ) 
i=i 9 i=i 81 

E 7 = u 7 + u 1 + ci 7 + d? + s 7 + s 7 + c 7 + c 7 + W + 6 7 , (2.15) 
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A 7 = —(3m 7 + 3m 7 + 3c 7 + 3c 7 - 2d? - 2<T - 2s 7 - 2s 7 - 26 7 - 26 7 ) . (2.16) 
15 

Because the photon is a charge conjugate eigenstate one can put the quark 
densities equal to the antiquark densities. 

We will now discuss the origin of the coefficient functions Ck,% (k = 2, L, 
i = q,g,j) which appear in (2.6). Starting with the last terms, the photonic 
coefficient functions Ck n are given up to next to leading order by the fol- 
lowing parton subprocesses. In the Born approximation the light quarks are 
produced by the reaction (fig. 2) 

T *(g) + j(k)^q + q, (2.17) 

while the heavy quarks are produced by the same reaction 

7*(q)+l{k)^Q + Q, (2.18) 

provided the square of the cm. energy denoted by s, where s — (k + q) 2 , 
satisfies the threshold condition s > 4m 2 . The 0(a s ) pQCD corrections are 
given by the one-loop contributions to processes (2.17) and (2.18) (see fig. 3) 
and the gluon bremsstrahlung processes (see fig.4) 

l*(q)+l(k)^q + q + g, (2.19) 

l*(q)+l(k)^Q + Q + g. (2.20) 
The parton cross section for the Born reaction in the case of light quarks 



2.17) can be found in fl4 |, [pSfl . In the case of heavy flavor production 



(2.18) the Born cross section is presented in [ IB] , p8| . Notice that the above 
reactions are very similar to the ones where the on-mass-shell photon 7(fc) is 
replaced by a gluon g{k). The cross sections of the photon-induced processes 
constitute the Abelian parts of the expressions obtained for the gluon-induced 
processes which are presented up to order a 2 s for the case of massless quarks in 



and in the case of massive quarks in |25|] . By equating some color factors 



equal to unity or zero in the latter expressions one automatically obtains 
the cross sections for the photon-induced processes above in particular for 
(2.19) and (2.20) (see Appendix). In the case of massless quarks the parton 
cross sections for (2.17), (2.19) contain collinear divergences which can be 
attributed to the initial photon being on-mass-shell. These singularities are 
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removed by mass factorization in the following way. We define 



f 1 f 1 Q 2 
Fk a ( z i Q 2 -. e ) = / dzi J Q dz 2 5(z - z 1 z 2 )T il (z 1 , M 2 ,e)C k)i (z 2 ,j^) , 

(2.21) 

where jF fc 7 (z, Q 2 , e) is the parton structure function, which is related to 
the parton cross section in the same way as the photon structure function 
F^(x,Q 2 ) is related to the cross section d 2 a/dxdy in (2.3). The parton 
structure function contains the collinear divergences represented by the pole 
terms e~i (j is a positive integer) where e = n — 4 (we use dimensional reg- 
ularization). These divergences are absorbed in the transition functions Tj 7 
(i — 7 , q , g) which depend both on e and on the mass factorization scale 



M. They can be inferred from the Abelian parts of T ig in [14|, [18], [29] and 



Both the photonic and hadronic coefficient functions Ck,i (i = 7 , q , g) 
which appear in the expressions for F^x^ 2 ) and F2(x,Q 2 ) in (2.6) are 
computed in the MS scheme. The coefficient functions C^k in (2.6) and 
(2.21) can be expanded in a power series in a s as follows 

ClM .cg» + 2^Ca» + (2fS) l Cg + ... (2.22) 

which holds for the light as well as the heavy flavor contributions. The pho- 
tonic coefficient functions for light quarks cf^ and can be directly derived 
via the mass factorization formula (2.21) from reactions (2.17) and (2.19) re- 
spectively. The heavy flavor coefficients and , which are obtained 
without using mass factorization, originate from processes (2.18) and (2.20). 
Notice that in the case of massive quarks the parton structure functions corre- 
sponding to the reactions (2.18) and (2.20) do not have collinear singularities 
and they can automatically be identified with the coefficient functions C^. 

Using the mass factorization formula in (2.21) one can also obtain the 
order a s contributions to the hadronic coefficient functions C^ q coming from 
process (2.19). The higher order contributions to the hadronic coefficient 
functions emerge when one calculates the NLO corrections to process (2.17). 
For example the gluonic coefficients Cj^ g can be inferred from the contribu- 
tions to 

l*(q)+l(k)^q + q + q + q } (2.23) 
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while g can be inferred from the contributions to 

7*(?)+7(*0 ^q + q + Q + Q- 



(2.24) 



Fortunately there is a quicker method to obtain the same information. The 
hadronic coefficient functions needed for the 0(a s ) renormalization group im- 
proved photon structure functions F^(x, Q 2 ) (2.6) can also be obtained from 
deep inelastic lepton hadron scattering, where the higher order corrections 
are known. For light flavor production we have listed the parton subprocesses 
and the corresponding coefficients which follow from these reactions in table 
1. We have given the corresponding information for heavy flavor production 
in table 2. In lowest order the photonic and hadronic coefficient functions 
have been presented in the literature (see [|14[ 



n, in 



J). Since these au- 
thors used a notation which is different from ours we will present the relevant 
formulae below. In next to leading order the expressions for the coefficient 
functions are obtained from (light quarks and gluons) and (heavy 
quarks) . However the expressions are too long to be presented in a paper Q. 
The method whereby the higher order coefficients can be derived from the 
expressions in ||24|| , |p5| is explained in the Appendix. 

Starting with the photonic coefficients for light quarks (see reaction (2.17)) 
they are given by 



Q_ 2 

M 2 ' 



Q_ 2 

M 2 



C%(z,^)=4{z 2 + (l-z) 2 }{ln^ + Ml 
and 



ln(z)} + 32z(l-z)-4, 
(2.25) 



Mo), 0^ 



1641 - z) . 



For massive quarks in the final state (see (2.18)) we have 



(2.26) 



C^°\z,Q 2 ,m 2 ) 



( 16m 2 , 

{4-8z(l -z) + ——z(l-3z) 



Q 2 



32m 4 



+ { - 4 + 322(1 - z) - 16^(1 - z))^l - — 



(2.27) 



1 These functions are available from smith@elsebeth.physics.sunysb.edu. 
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and 



C»f(z,Q 2 ,m 2 ) = 16z(l-z) 



4m 2 m 2 

— A ±j 

s s 



(2.28) 



where m is the heavy-flavor mass and ^fs is the cm. energy of the virtual 
photon-real photon system. Furthermore we have 



'X-z 



Q 2 



L = ln 



1 + Jl -Am 2 /s 



4m 2 / s 



(2.29) 



Formulae (2.27) and (2.28) can be found in 

In the next order in a s process (2.19) (Fig.4) and the one-loop correc- 
tions to process (2.17) (Fig. 3) give rise to the coefficients C ( ^(z , Q 2 / M 2 ) . 
In the case the outgoing fermion lines in figs. 3,4 stand for the heavy flavors 
(see reactions (2.18) and (2.20)) the corresponding coefficients are given by 
Cj?!y (z, Q 2 , m 2 ). More information about the higher order corrections to the 
photonic coefficient functions can be found in the Appendix. 

In zeroth order of a s the hadronic coefficient functions are 



M0) ( Q 2 

C {0) (z ^~ 
M0) ( 0^ 



5(1 

0, 
0. 



(k = 2,L). 



(2.30) 
(2.31) 
(2.32) 



In order as the hadronic coefficient functions originating from a light quark 
in the initial state (table 1) are given by 



M 2 ' 



C h 



-2-2z} 



x{l « +M l-,)-2} 



M 2 



n 1 + z 2 , 9 5 

2 In 2 H V-z 

l-z 2 2 



-5(l-z){3\n 



M 2 



-9 



4C(2)} 



where the standard definition of a plus distribution is used, and 

Q 2 



c£l(z, 



M 2 



43 



(2.33) 



(2.34) 
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Notice that in order a s there is no difference between Cf'^ and C^g . The 
coefficient functions for a gluon in the initial state and massless quarks in the 
final state (table 1) can be derived from (2.25) and (2.26) via multiplication 
by a color factor 



4%,j^)-n f T f C^(z^), (k = 2,L). (2.35) 



An analogous relation holds when the massless quarks in the final state are 
replaced by the heavy flavors (table 2) and we get from (2.27) and (2.28) 

C^ 1 \z,Q 2 ,m 2 )=T f C k °! / (z,Q 2 ,m 2 ), (k = 2, L) . (2.36) 

The color factors which appear in the above equations are given by Cf = 4/3 
and Tf = 1/2 for the case of SU(3). 

The higher order a 2 corrections to the coefficient functions, describing 
massless partons only are denoted by C ki where i = q,g (see table 1). They 



have been calculated in 24 . In the Appendix we have decomposed Cl [ into 



color factors so that we can infer the 0(a s ) photonic coefficients C k ^ from 

the Abelian part of C kg . 

The 0(a 2 ) corrections to the heavy flavor coefficient functions given by 
C ( k J(z,Q 2 /M 2 ,m 2 ) and C k f\z, Q 2 /M 2 , m 2 ) (table 2) are calculated for the 
first time in |25l . The relations between these coefficients and the ones derived 



in section 5 of will be presented in the Appendix. By decomposing them 
in color factors we again can derive the photonic heavy flavor coefficient 
Cjfi from the Abelian part of C kg . Since in lowest order the hadronic 

heavy flavor coefficient C^ k 2 }(z,Q 2 /M 2 ,m 2 ) only contributes up to the 0{a 2 ) 
level, when i = q we do not have to distinguish between singlet (S) and 
non-singlet (NS) and we can put 



C S k f(z, ^m 2 ) = C (2) (^ ^rn 2 ) = C%(z,Q 2 ,m 2 ) . (2.37) 



The above expression indicates that in lowest order C k q (z,Q 2 ,m 2 ) is deter- 
mined without having performed mass factorization which is indicated by its 
independence of the mass factorization scale M. This is because it originates 
from the Compton scattering process, which in lowest order does not have 
collinear singularities. 
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Finally in table 3 we have translated our notations for the coefficient 
functions into those used in [|14]], [p3|| , |28|| , fp9(1 . We also list the new contri- 



butions to the photon structure functions which were not included earlier in 
the literature. 
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3 Results 



In this section we will discuss the NLO QCD corrections to the photon struc- 
ture functions (x, Q 2 ) for k = 2, L. In particular we focus our attention 
on the heavy flavor contributions (mainly charm), which originate from the 
hadronic as well as the photonic coefficient functions in (2.6). Since heavy 
flavors can be produced either in virtual-photon parton or in virtual-photon 
real-photon reactions we will call the former hadronic heavy flavor production 
and the latter photonic heavy flavor production. 

In the subsequent part of this section we want to make a comparison 
between the LO and NLO description of the photon structure functions, 
where all contributions listed in tables 4 and 5 are included. Furthermore we 
want to investigate the relative magnitude of the heavy flavor (mainly charm) 
component of the structure function. We also show the difference between 
the massless and massive heavy flavor approach. When the heavy quarks are 
treated as massless, their contribution to the photon structure functions are 
given by the corresponding parton densities in the photon convoluted with the 
light quark and gluon coefficient functions. This description is appropriate 
when Q 2 » m 2 . If Q 2 is of the same order of magnitude as m 2 , then the 
massive quark approach has to be adopted and the heavy flavor production 
is described by the heavy flavor coefficient functions in (2.6) which can be 
computed order-by-order in perturbation theory. 



In the literature a LO analysis was given for F%{x, Q ) in [16] and F2(x, Q 



m 



23fl . Here all LO coefficient functions in tables 4 and 5 were included ex- 



cept for the ones related to hadronic heavy flavor production, (i.e., 7* +g — ► 
Q + Q). The last contributions were also neglected in the NLO analy- 
sis for F2(x,Q 2 ) in pi| and the photonic heavy flavor contribution from 
7* + 7 — > Q + Q was only taken into account in lowest order. A NLO anal- 
ysis of F2(x,Q 2 ) could not be carried out previously because the order a 2 
contributions to all the longitudinal coefficient functions were not known un- 
til recently. Since all NLO coefficient functions are now known, and they are 
listed in tables 4 and 5, we are able to present a complete NLO description 
for both F^{x, Q 2 ) and for F2(x, Q 2 ) as well as make a comparison with the 
LO descriptions. 

In our plots we adopt the LO and NLO parametrizations of the parton 
densities in the photon from |19j (for other sets see | 20| , pl|). For rif = 



3 we use Aqcd = 232 MeV at leading order and Aqcd = 248 MeV at 
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next to leading order. For n/ = 4, both the leading order and the next 
to leading order Aqcd are set equal to 200 MeV. In leading order, we use 
a one-loop result for the running coupling constant and in next to leading 



order a two-loop corrected running coupling constant is chosen, see e.g. |L£ 
All calculations are done with M 2 = Q 2 , except where otherwise indicated. 
In our analysis, when we take the charm quark to be massive, we take m c = 
1.5GeV/c 2 . Furthermore we take three light flavors (rif — 3) for the parton 
densities, the coefficient functions and the running coupling constant. When 
we treat the charmed quark as massless, it then takes on the identity of an 
ordinary parton, so we set rif = 4. The bottom and top quark contributions 
will be omitted since they are negligible for the Q 2 values accessible at past 
and present experiments. In the LO approximation the corresponding parton 
densities are multiplied by the coefficient functions in tables 4 and 5, which 
are indicated by LO. In NLO we have chosen the MS scheme for the parton 
densities, the coefficient functions and the running coupling constant. The 
coefficient functions which have to be added to the LO ones are indicated 
by NLO in tables 4 and 5. In order to get a consistent NLO analysis for 
the structure functions we follow the procedure in |l||, which is explained 
in Therefore we multiply the LO coefficient functions by p and the 

NLO coefficient functions by p in (2.6) (for the notation of p and p see 
eqn. (A. 23) and the discussion in the Appendix A in |19||). Notice that in 



IS] the parton densities described in Appendix A were presented in the DIS„ 



scheme. However they can be changed into the MS scheme via eqns.(4)-(6) 
in [19]. After changing the lowest order photonic coefficient funct ion C^/y in 
the DIS 7 scheme we have checked that both schemes lead to the same result 
provided the change of eqn. (4) in [18| is only applied to the parton density 
denoted by p as defined above. 

We now compare the results from our calculations for F^(x, Q 2 ) first with 
data from PLUTO Q (Q 2 = 5.9 (GeV/c) 2 ) and then with data from AMY 
(Q 2 — 51 (GeV/c) 2 ). We also show predictions for F2{x,Q 2 ). 

In fig. 5 we make a comparison between the LO and NLO aproximation 
for i^x, Q 2 ) at Q 2 = 5.9(GeV/c) 2 , where the heavy charm components 
(hadronic and photonic) are included. The low-x hump is due to charm 
production, which turns off at about x = 0.4 (the threshold value). We also 
show separately the contributions due to massive charm production. When 
this contribution reaches its maximum value it constitutes about 20 % of 
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the structure function F^ in LO and 30 % in NLO. The 0(as) correction to 
the Born contributions to massive charm production are quite large, adding 
approximately another 50 % to the Born terms. Overall, we observe that LO 
and NLO are not very different. Note that the data also seem to indicate the 
presence of a charm component. 

In fig. 6 we do the same for F2(x,Q 2 ) at Q 2 = 5.9 (GeV/c) 2 . This is for 
theoretical purposes only: there are no data presently available for F2 at any 
value of Q 2 . We see from this and the previous figure that there is not much 
difference between the LO and NLO results both for K7 and F2- However, 
the heavy charm component of F2 is less important than in the case of K?. 
At LO it is about 15 % where this component reaches its maximum, whereas 
in NLO it amounts to about about 30 % also. The latter is due to the fact 
that the 0(a s ) corrections to the heavy charm component of F2 are as large 
as 100 %. 

In fig. 7 we present F2(x,Q 2 ) at LO for three different choices of mass 
factorization scale. Note that in this case the only variation is due to the 
parton densities. The variation in the M dependence is uniform over the 
whole x-range. In fig.8 we do the same for F2(x, Q 2 ) at Q 2 = 5.9 (GeV/c) 2 . 
Here there is additional scale dependence due to a s (M 2 ). Hence, contrary 
to fig. 7, the curve for M = Q/2 is the upper one here. 

Fig. 9 shows the same as fig. 7 but now at NLO. There is now additional 
scale dependence due to a s (M 2 ) and the mass factorization scale logarithms 
of the type \n(Q 2 /M 2 ) in the coefficient functions (see e.g (2.25) and (2.33)). 
Note that the scale dependence is reduced in the small-x region compared to 
the LO case. However at very large x values, where the charm contribution 
can be neglected, the scale variation is larger than in the LO case. This 
is due to the pointlike light quark contribution, which drops increasingly 
dramatically as one increases M. At small x this is partially offset by the 
increase of the charm contribution. 

In fig. 10 we show the same plots as in fig. 9 for F2(x,Q 2 ). The scale 
variation is small as in the LO case. 

We now turn to a comparison of results for massive versus massless charm 
contributions as defined above. Since the differences are essentially the same 
in the LO case as in the NLO case we only show plots for the latter. Therefore 
in fig. 11 we compare the NLO massless (rif = 3) plus the massive charm- 
quark contribution to F^x^Q 2 ) at Q 2 = 5.9 (GeV/c) 2 , with the NLO mass- 
less (rif = 4) contribution. Note that the massless rif — 4 contribution is 
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smaller than the curve where we take rif = 3 massless and a massive charm 
quark, even at large x where the charm contribution is zero. This is due 
to a change in Aq C d and consequently a change in the parton distribution 
functions. However the difference between the massless and massive cases is 
small for large x, where threshold effects are negligible. 

In fig. 12 we show the same plots for F2(x, Q 2 ) in NLO at Q 2 = 5.9 (GeV/c) 2 . 
Note the enormous increase that occurs in going to the rif — 4 massless case. 
Since this effect is already there in the LO case it can be understood as fol- 
lows. In the case of rif = 3 where charm is considered massive one includes 
the coefficient function (2.26), which is multiplied by 2/9 and 
(2.28), which is multiplied by 16/81. If rif — 4 the charm is treated as mass- 
less and C 1 ^'^ (massive charm) is replaced by (massless charm). Since 
the latter is much larger than the former due to the additional suppression 
factor in (2.28) this explains why the result for rif = 4 is much larger than 
for rif = 3. 

In fig. 13 we show the ^-dependences of the massive hadronic charm con- 
tribution and the massive photonic charm contribution to F^(x, Q 2 ) at Q 2 = 
5.9 (GeV/c) 2 in LO and in NLO. The corresponding results for F2(x, Q 2 ) are 
shown in fig. 14 in LO and in NLO. The interesting feature to note in all these 
figures is the complete dominance of the photonic charm production over the 
hadronic production. This makes F2{x,Q 2 ) for massive charm production 
at moderate x a very promising test of pQCD, because of the lack of depen- 
dence on the hadronic component. Experimentally this is of course a very 
difficult quantity to determine, but perhaps not impossible. The same holds 
for F2(x, Q 2 ) for massive charm production, but that is even more difficult to 
determine experimentally. However for x < 0.01 the pointlike contributions 
to both F~2 and F2 for massive charm production become very small and the 
hadronic component begins to dominate. 

We now repeat all the figures for the Q 2 = 51 (GeV/c) 2 value of the AMY 
collaboration. We remark that now the charm contribution switches off at 
x = 0.85. Here the heavy charm component becomes in general larger than 
in the case for Q 2 = 5.9(GeV/c) 2 . For F 2 7 it is 30 % in LO where this 
component reaches its maximum, and 40 % in NLO. For F^ the percentages 
are roughly similar. Note however that the 0(a s ) corrections are smaller 
than for Q 2 = 5.9 (GeV/c) 2 . For F 2 7 they are up to 15 % and for Fl up 
to 30 %. The mass factorization scale dependence at large x for F^{x, Q 2 ) 
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at NLO seems (fig. 19) to be somewhat reduced compared to the case of 
Q 2 = 5.9 (GeV/c) 2 but still larger than at LO (fig. 17). 

To conclude, we have presented in this paper the first complete NLO 
analysis of F^{x, Q 2 ) and F2{x, Q 2 ) containing both light and heavy quarks. 
Summarizing our findings we have seen that for both values of Q 2 we con- 
sidered the NLO structure functions are not too different from the LO ones. 
This is not so surprising for F^ (x, Q 2 ) since we used the parton densities of 
19|j and most of the contributions were already included in their analysis ex- 



cept for O(as) corrections to heavy quark production, which are numerically 
small. We see that F% has a moderate sensitivity to changes in the mass 
factorization scale except at large x. 

For F^{x, Q 2 ) this is the first NLO analysis, and at the same time com- 
plete, since all heavy and light quark contributions have been included. We 
found that F^{x, Q 2 ) changes very little from LO to NLO, and is very stable 
under scale changes. Above x ~ 0.1 the hadronic production of charm is 
small compared with the photonic production, while the former is dominant 
for x < 0.01. All this would make a measurement of F2(x,Q 2 ) (e.g. at 
LEP2) an interesting prospect. 

Our results could be used to determine more accurate NLL parton dis- 
tribution functions for the photon. This would become especially relevant 
when data become available for for charm production, and for F2- Fi- 
nally, we stress that if the heavy quark contribution could be extracted from 
a measurement of F% this would yield a very good test of perturbative QCD. 
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Table 1. 



order 


parton subprocess 


coefficient function 


«2 


7* + q(q) -»■ 5(5) 


r (0) 




7* + g(g) -»• 9(9) + 5 


r NS,(l) r s,(i) 

^k,q - ^k,q 




7* + -»• g + g 


C k,g 


a 2 


7* + 5(5) -»■ 5(5) + 9 + 9 


r NS,(2) r S,(2) 
^k,q - ^k,q 




7* + q(q) -»• 9(9) + 9(9) + 9(9) 


r NS,(2) , p S,(2) 
u fc,g T ^k,q 




7* + #^9 + 9 + 5' 


V k,g 



List of deep inelastic virtual-photon-parton subprocesses up to 0(a 2 ). 
The one and two-loop corrections to the lower order processes have been 
included in our calculations but are not explicitly mentioned in the table. 

Table 2. 



order 


parton subprocess 


coefficient function 


«! 


7* + 9 -> Q + Q 




a 2 s 


l* + 9^Q + Q + 9 


U k,g 


a 2 


7* + 9(9) -> 9(9) + Q + Q 


n H,(2) r NS,(2) r S,(2) 
°fc,g > k,q — ^k,q 
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List of deep inelastic virtual-photon-partonic subprocesses contributing to 
heavy flavour production up to 0(a 2 s ). The one-loop corrections to the Born 
approximation have been included in our calculations but are not explicitly 
mentioned in the table. 

Table 3. 



this paper 




M 


23 * , 28 ** 


(2.25) 


(4.12) 


5 7 (3.7) 




(2.26) 






(15)* 


^(§) 4 C£f (2.27) 






\Flc (2-13)** 


iKfXf (2.28) 








C% (2-33) 


(4-10) 


Bm,B q (3.7) 




C£i (2-34) 








C?J (2.35) 


4 nJ (4.11) 


5 G (3.7) 




C (2-35) 








C lj ( 2 - 36 ) 










24 1 












2r 


1 








cSS I- 


24 1 












2£ 
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Notations in several papers for the hadronic and photonic coefficient func- 
tions. Notice that the expressions in [2S] are in Mellin transform space. 
The blanks mean that these contributions were not considered in the papers 
quoted. 

Table 4. 



order 


parton subprocess 


coefficient function 


a 


7* + <?(<?) -> <?(<?) 


CgJ LO, NLO 




7* + 7 — > q + q 


4°} NLO 






C 2 ^ 0) LO, NLO 


a] 


7* + q(q) -> q(q) + g 


C w (=c?i 1J ) NL0 




7* + g -> g + g 


Cg NLO 




7* + -> Q + 


C£/ 1} LO, NLO 




7* + 7^Q + Q + 


c£; w NLO 




7* + 0^Q + Q + 


NLO 




7* + q(q) -> + Q + 


C a, ,C ia, (=0 NL0 



Coefficient functions used in this paper for a leading order (LO) and a 
next-to-leading order (NLO) analysis of F%(x, Q 2 )/a. 

Table 5. 
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order 


parton subprocess 


coefficient function 




7* + 7 — > g + g 


LO, NLO 




7*+7~*Q + Q 


C"f } LO, NLO 




7* + -> + g 


C ,w (=0 L0 > NL0 




7* + # -> q + q 


C?] LO,NLO 




i* + g^Q + Q 


C?j 1J LO, NLO 




Y + i^q + q + g 


NLO 




7* + 7 ^Q + Q + g 


NLO 


a 2 


7* + + # + # 


C' W (=0 NLO 




7* + 9(9) q(q) + q(q) + q(q) 


C^CVC") NLO 




Y + g^q + q + g 


NLO 




T + g-^Q + Q + g 


Cf i9 (2) NLO 




7* + q(q) -> + Q + Q 





Coefficient functions used in this paper for a leading order (LO) and a 
next-to-leading order (NLO) analysis of F2(x,Q 2 )/a. 



Appendix 

In this Appendix we show how one can derive the 0(a 2 ) coefficients corre- 
sponding to the reactions in tables 1 and 2 from the expressions calculated 
in j2J| and [23] respectively. The 0(a 2 ) coefficients mentioned in table 1 are 
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given by 



NS,{2) ( Q 2 v _ r 2nW/ Q 2 N , ^ r H (*)/_ ^ 

*,« ^' ~~ ^F^ffKZ, j^) + u a ^f^af{ z , Jp. 



+ % T F ^4^(z, ^) , (A.l) 

where andC^'^ are the coefficients of the (« s /47r) 2 term in eqns.(B.l) 

and (B.2) of |24| respectively. The singlet coefficients can be split into a non- 
singlet and a pure singlet piece as follows 

The pure singlet coefficients Cjf?'^ can be written as 



C (2) (*> %) = n f T f C F DP F (z, ^) , (A.3) 



where C^ S '^ and C^q are the coefficients of the (a s /4ir) 2 terms in eqns.(B.3) 
and (B.4) of |24| respectively. Finally the gluonic coefficient is given by 

0*> Jjs) = n f T f C F Ef F {z, + n f T f C A Ej»(z, ^) , (A.4) 

where C^ g and C| 2 g are the coefficients of the (a s /47r) 2 terms in eqns.(B.5) and 
(B.6) of [|2|] respectively. The color factors in SU(3) are given by Cp = 4/3, 
Ca — 3, Tp = 1/2 and riy denotes the number of light flavors. The 0(a s ) 
photonic coefficient can be derived from the Abelian part of C^ g (A.4) 
and it equals 



O*. ^) = ^) ■ (A.5) 



The coefficient functions due to heavy flavor production (see table 2) are 
related to the coefficients defined in |25] in the following way. In first order 
in a s we have (see also (2.36)) 

C?fWV) = -^ci>,0, (A.6) 
25 



with 



C?f\z,Q 2 ,m 2 ) = ^{cS^O + ^O}, (A.7) 

/I fit z, 



»? = t £ i-i > e = ^- (a.8) 

4m z rrr 



In second order in a s one gets for i = q,g 



Cg> (*, Q 2 , m 2 ) = l&r-^-d&fo, , (A.9) 



Cg(z,Q 2 ,m 2 ) = 167T J^{4>,£) H-dg^O}, (A.10) 



and 

CS (2) (*, 7H>™ 2 ) = 16 ^i c S(^0 + €(V,0^} , (A.ll) 



2 



C?f\z, gr, m 2 ) = I67r^{cgfo, + ^(v, + Ngfo, 

+#,ifo,0]ln^}. (A.12) 
rrr 

In the above expressions the coefficients c^ 1 ], cjj. 1 ] and g^ 1 ] for k = T,L and 
i — q,g are defined in eqns.(5.3)- (5.6) of |[25|| . As has already been mentioned 
they are too long to be presented in a paper and they are available upon 
request. Like the coefficient functions in table 1 the heavy flavor contributions 
can be decomposed in color factors in a similar way. In first order in a s we 
have 

Cff\z,Q 2 ,m 2 ) = T f Cff(z,Q 2 ,m 2 ) , (A.13) 

where C^'^ denotes the photonic coefficient which is given in eqs.(2.27) and 
(2.28) (see also (2.36)). In second order in a s the expressions are analogous 
to the ones presented for light quark production in (A.l), (A, 3) and (A. 4) 

ci%,Q 2 ,m 2 ) =T f C F BP F (z,Q 2 ,m 2 ), (A.14) 
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C 2) (*> ^V) = T/C7 F £>S (z, ^,m 2 ) , (A.15) 

and 

^V) = TfCpE^l(z, ^,m 2 ) + T^i^, |^,m 2 ) .(A.16) 

Notice that in the limit m — > the above expressions need an additional 
mass factorization. After this procedure is carried out the coefficients Bpp, 
Dff, Eff and Efa pass into their massless analogues defined in (A.l), (A. 3) 
and (A. 4). The order a s contributions to the photonic coefficient function 
can be derived from (A.16). It is equal to 

C^\z,Q\m 2 ) = C F EP F (z,Q 2 ,m 2 ), (A.17) 
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Figure Captions 

Fig.l. The process e~{p e ) + e + — > e~(p' e ) + e + + X, where X denotes any 
hadronic state. 

Fig. 2. The lowest order Feynman diagrams contributing to the Born reac- 
tion 7*(g) + 7(fc) — > q + q. 

Fig. 3. Feynman diagrams contributing to the one-loop correction to the 
process 7*(?) + 7(fc) q + q. Additional graphs are obtained by 
reversing the arrows on the quark lines. Graphs containing the external 
quark self-energies are included in the calculation but not shown in the 
figure. 
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Fig. 4. The order g (a s = g 2 /47r) Feynman diagrams contributing to the 
gluon bremsstrahlung process 7* (5) + j(k) — > q + q + g. Additional 
graphs are obtained by reversing the arrows on the quark lines. 

Fig. 5. The x-dependence of F^(x,Q 2 ) at Q 2 = 5.9 (GeV/c) 2 , solid line: 
F^iNLO), long-dashed line: i^(LO), short-dashed line: NLO heavy 
quark contributions, dotted line: LO heavy quark contributions. The 
data are from PLUTO 0. 

Fig.6. The x-dependence of F2(x,Q 2 ) at Q 2 = 5.9 (GeV/c) 2 , solid line: 
F2{NLO), long-dashed line: F2(LO), short-dashed line: NLO heavy 
quark contributions, dotted line: LO heavy quark contributions. 

Fig.7. The x-dependence at LO of F 2 7 (x, Q 2 ) at Q 2 = 5.9 (GeV/c) 2 for three 
choices of the mass factorization scale M 2 : M = 2Q (long-dashed line), 
M = Q (solid line) and M = Q/2 (short-dashed line). The data are 
from PLUTO 0. 

Fig.8. The x-dependence at LO of F^(x, Q 2 ) at Q 2 = 5.9 (GeV/c) 2 for three 
choices of the mass factorization scale M 2 : M = 2Q (long-dashed line), 
M = Q (solid line) and M = Q/2 (short-dashed line). 

Fig.9. The x-dependence at NLO of F 2 7 (x, Q 2 ) at Q 2 = 5.9 (GeV/c) 2 for 
three choices of the mass factorization scale M 2 : M = 2Q (long-dashed 
line), M = Q (solid line) and M = Q/2 (short-dashed line). The data 
are from PLUTO |]. 

Fig.10. The x-dependence at NLO of F 7 (x,Q 2 ) at Q 2 = 5.9 (GeV/c) 2 for 
three choices of the mass factorization scale M 2 : M = 2Q (long-dashed 
line), M = Q (solid line) and M = Q/2 (short-dashed line). 

Fig. 11 The x-dependence of the NLO massless (n/ = 3) plus the mas- 
sive charm-quark contribution to F^(x, Q 2 ) (solid line) compared with 
the NLO massless contribution (n/ = 4, dashed line), at Q 2 = 5.9 
(GeV/c) 2 . The data are from PLUTO 0. 

Fig. 12. The x-dependence of the NLO massless (n/ = 3) plus the mas- 
sive charm-quark contribution to F^{x, Q 2 ) (solid line) compared with 
the NLO massless contribution (n/ = 4, dashed line), at Q 2 = 5.9 
(GeV/c) 2 . 
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Fig. 13. The x- dependence of the LO and NLO massive hadronic charm 
contributions to F^ (x, Q 2 ) (solid lines) compared with the LO and 
NLO massive photonic charm contributions (dashed lines), at Q 2 = 5.9 
(GeV/c) 2 . The NLO contributions are the larger ones. 

Fig. 14 The x-dependence of the LO and NLO massive hadronic charm con- 
tributions to F2(x, Q 2 ) (solid lines) compared with the LO and NLO 
massive photonic charm contributions (dashed lines), at Q 2 = 5.9 
(GeV/c) 2 . The NLO contributions are the larger ones. 

Fig. 15. The x-dependence of K^x, Q 2 ) at Q 2 = 51 (GeV/c) 2 , solid line: 
F2(NLO), long-dashed line: F2(LO), short-dashed line: NLO heavy 
quark contributions, dotted line: LO heavy quark contributions. The 
data are from AMY || . 

Fig.16. The x-dependence of F2(x,Q 2 ) at Q 2 = 51 (GeV/c) 2 , solid line: 
F2(NLO), long-dashed line: F2(LO), short-dashed line: NLO heavy 
quark contributions, dotted line: LO heavy quark contributions. 

Fig.17. The x-dependence at LO of F 2 7 (x, Q 2 ) at Q 2 = 51 (GeV/c) 2 for three 
choices of the mass factorization scale M 2 : M = 2Q (long-dashed line), 
M = Q (solid line) and M = Q/2 (short-dashed line). The data are 
from AMY [§. 

Fig.18. The x-dependence at LO of F2(x, Q 2 ) at Q 2 = 51 (GeV/c) 2 for three 
choices of the mass factorization scale M 2 : M = 2Q (long-dashed line), 
M = Q (solid line) and M = Q/2 (short-dashed line). 

Fig.19. The x-dependence at NLO of F 2 7 (x,Q 2 ) at Q 2 = 51 (GeV/c) 2 for 
three choices of the mass factorization scale M 2 : M = 2Q (long-dashed 
line), M = Q (solid line) and M = Q/2 (short-dashed line). The data 
are from AMY 0. 

Fig.20. The x-dependence at NLO of F2(x,Q 2 ) at Q 2 = 51 (GeV/c) 2 for 
three choices of the mass factorization scale M 2 : M = 2Q (long-dashed 
line), M = Q (solid line) and M = Q/2 (short-dashed line). 

Fig. 21 The x-dependence of the NLO massless (n/ = 3) plus the massive 
charm-quark contribution to F% (x, Q 2 ) (solid line) compared with the 
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NLO massless contribution {rif = 4, dashed line), at Q 2 = 51 (GeV/c) 2 . 
The data are from AMY || . 

Fig. 22. The x-dependence of the NLO massless (rif = 3) plus the massive 
charm-quark contribution to F2{x,Q 2 ) (solid line) compared with the 
NLO massless contribution {rif = 4, dashed line), at Q 2 = 51 (GeV/c) 2 . 

Fig. 23. The x-dependence of the LO and NLO massive hadronic charm 
contributions to (x, Q 2 ) (solid lines) compared with the LO and 
NLO massive photonic charm contributions (dashed lines), at Q 2 = 51 
(GeV/c) 2 . The NLO contributions are the larger ones. 

Fig. 24 The x-dependence of the LO and NLO massive hadronic charm con- 
tributions to F2(x,Q 2 ) (solid line) compared with the LO and NLO 
massive photonic charm contributions (dashed lines), at Q 2 = 51 (GeV/c) 2 . 
The NLO contributions are the larger ones. 
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